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Abstract. Searches for lepton flavor violating r l-y, r — > U£ and r — > Ihh decays 
at the BABAR experiment are presented. Upper limits on the branching ratios are 
obtained at the level of O(10~^) at 90% confidence level. 

1 Introduction 

Lepton flavor violating (LFV) processes such as the neutrinoless decay of the 
T lepton have long been identified as unambiguous signatures of new physics, 
because no known fundamental local gauge symmetry forbids such a decay. 
While forbidden in the Standard Model (SM) because of vanishing neutrino 
mass in the three lepton generations, extensions to include current knowledge 
of neutrino mass and mixing imply B{t —^ ^j) ~ 0(10^^^) [1], which is many 
orders of magnitude below the experimental sensitivity. However, many new 
theories, as tabulated below, allow for LFV decays: r £7, t U£, t Ihh 
(where £ = e, ^; h = ir, K) up to their existing experimental bounds ^ 0(10^^): 
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Feynman diagrams for r /i7 and t — > /i/i// decays via s-neutrino mixing in 
minimal supergravity model with heavy z^r (seesaw mechanism) and via neutral 
Higgs exchange in supersymmetric model are shown in Figure 1, respectively. 





Figure 1: Illustrative scenarios for r — > /i7 (left) and r —t /i/i/i (right). 
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2 Signal Identification 



Searches for LFV decay in modes: r ^ ^7 [1,10], r Ui [11] and r ihh [12] 
have been performed with 232.2 fb~^, 91.6 fb~^ and 221.4 fb^^ of data collected 
by the BaBAR experiment at y/s « 10.58 GeV, respectively. The characteristic 
feature of these decays is that both the energy and the mass of the r-daughters 
are known in such an e^e~ annihilation environment. In terms of the two 
independent variables: beam energy constrained mass (mEc) and the energy 
variable AE = Et- — y/s/2, where Et is the energy of the r-daughters in center- 
of-mass system, the signal is clustered around (mr,0) in the (mEc, AE) plane. 

The identification of daughters from signal r decays are optimized for searches 
for each decay mode separately. The electrons are identified from the energy 
deposited in the electromagnetic calorimeter and momentum of the track mea- 
sured in the drift chamber with an efficiency of 91% for r — > 57 and r ££i 
searches. The muons are identified by its minimal ionizing particle signature 
in the calorimeter and hits in the instrumented fiux return with an efficiency 
of 82%, 63% and 44% in r fij, r £££ and r — > £hh searches respectively. 
The kaons are identified using the measured rate of ionization loss in the drift 
chamber and the measured Cherenkov angle in a ring-imaging detector with an 
efficiency of 81%. The mis-identification rates for a pion track to be identified 
as an electron, a muon or a kaon are 0.1%, 1.0 — 4.8%, 1.4% respectively. 

3 Background estimation 

The primary backgrounds are from Bhabha or di-muon events, which are re- 
stricted to a narrow band at small values of |AE|, or the e+e^ — > r+r^ events, 
which are restricted to negative values of AE, because the signal topology 
reconstruction does not account for the missing neutrino's. The remaining 
backgrounds from e+e^ — )■ qq are uniformly distributed. 




Figure 2: mEC distribution inside a ±2(t band in AE for t ^ £j searches. 



For T ^ £j searches, the signal probabihty density function (PDF) is de- 
scribed by a double Gaussian shape in niEc, and the background is well de- 
scribed by a constant PDF or with a small slope in niEc inside a ±2cr band in 
AE, as shown in Figure 2 for t ej (left) and t ^-f (right) decays. 

For T III (Ihh) searches, the background PDF's are analytically parame- 
terized as function of AE and AM = uidau - m,- , where mdau is the reconstructed 
mass of the t daughters. The background rates are determined by un-binned 
maximum likelihood fits to the data, shown along with the selected signal MC 
events in Figure 3 and 4 for t — > III and r Ihh decay modes respectively. 

All these searches are performed in a blinded manner, where the background 
predictions from sideband data are compared to the data inside the signal 
region, only after the optimization and systematic studies have been completed. 
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Figure 3: Observed data as dots and the boundaries of the signal region in 
(AM, AE) plane for t III searches. The dark and light shading indicates 
contours containing 50% and 90% of the selected MC signal events, respectively. 
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Figure 4: Observed data as dots and the boundaries of the signal region in 
(AM, AE) plane for r —^ Ihh searches. The dark and light shading indicates 
contours containing 50% and 90% of the selected MC signal events, respectively. 



4 Results 



No signal has been observed. Upper limits at 90% confidence level (C.L.) are 
set using: Bf^^^ = Nfj\/ {2£C(Jrr), where Nf^^ is the 90% C.L. upper hmit on 

the number of signal events for iVobs events observed when A^bgd background 
events are expected, and e is the signal efficiency. Efficiency estimates, the 
number of expected background events (A^bgd) in the signal region (with total 
uncertainties), the number of observed events (A'obs) in the signal region, and 
the 90% C.L. uppcT limit {Bf^-i) for each dcx;ay modc^ arc> tabulatc^d below: 
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5 Summary 

An improvement of five order of magnitude in the upper limits for 4 LFV r 
decay modes over the past two decades is shown in Figure 5. The next five years 
promises to be most interesting phase in this evolution, when experimental 
results approach closer the predictions from different theoretical models. 
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Figure 5: Evolution of experimental bounds (Sul) ^-nd some predictions. 
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